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ABSTRACT 

We present results on the clustering of Lyman break galaxies (LBGs) around a probable damped Lya 
absorption (DLA) line cloud at Zabs — 2.974 from deep UBVI images of the field containing the quasar 
APM 08279+5255 (zcm — 3.91). The large area covered by our images, 0.31 deg^ or ~ 40 X 40 Mpc 
co-moving at redshift z = 3, and their depth, ni^xaisky) ~ 27.6 mag arcsec"^, allow us to identify 
~ 450 LBG candidates brighter than Iab = 24.80'at 2.75 < z < 3.25 both close (50 kpc) to the DLA 
line-of-sight and up to 20 Mpc (co-moving) from the DLA, i.e. physically unrelated. LBG candidates 
were identified using photometric redshift techniques that include the / magnitude as a prior estimate in 
addition to the colors. The two are combined using Bayes' theorem. This helps to break the degenerancies 
that occur in a pure spectral template fitting scheme. The overall rms is (Tz ~ 0.15 at z 3 based on 
our analysis of photometric redshifts in the HDF-N. 

From the redshift likelihood distributions, we selected LBG galaxies within a redshift slice of width 
Wz = 0.15(~ az) centered on the redshift of the DLA Zabs- 

Within that redshift slice, we find an enhancement of galaxies near the DLA using both the surface 
density (S/Sg ~ 3) and an estimator of the 3-D spatial overdensity {n/rig ^ 5 ± 3). The surface 
overdensity S/Sg is significant at the > 95% significance level on scales 2.5 < rg < 5 Mpc co-moving. 
The overdensity can not be related to the QSO environment since the QSO is at Zem = 3.91. These 
results imply that some DLA could reside in high density regions. 

We search within 45" from the line of sight for galaxies responsible for the DLA, and find one candidate 
with Zphot = 3.03 that is 26" (145 kpc physical) away. From its magnitude / — 24.65 + 0.2, its luminosity 
is M/^AB — —21.35. Due to its large impact parameter, however, this galaxy is not a likely candidate for 
the absorber. 

Subject headings: galaxies: evolution — galaxies: high-redshift — quasars: absorption lines — quasars: 
individual (APM08279+5255) 
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1. INTRODUCTION 



Damped Lya absorbers (DLAs) give rise to the largest 
neutral hydrogen (Hi) column densities (Ajjj > 2 x 
10^° cm~^) of quasar absorption line systems. They are 
important because they contain the largest reservoir of Hi 
at high redshift, and their detailed nature have been the 
subject of an ongoing debate. Wolfe et al. (1986) argued 
that DLAs are massive gas rich disks based on the similar- 
ity between the amount of Hi contained in DLAs and the 
stellar masses of disk galaxies today. Later, Wolfe et al. 
(1995) and Prochaska & Wolfe (1997) argued that DLA 
kinematics are consistent with those expected from lines 
of sight intercepting massive rotating thick gaseous disks. 

However, theoretical simulations of galaxy formation 
show that a large range of structures and morphologies can 
give rise to DLAs (e.g. Katz et al. 1996; Haehnelt, Stein- 
metz & Ranch 1998; McDonald & Miralda-Escude 1999). 
DLAs would arise from the combined effect of a massive 
central galaxy and/or a number of smaller satellites in 
a virialized halo (Mailer et al. 2000) or filaments (e.g. 
Haehnelt, Steinmetz & Ranch 1998), i.e., DLAs would 
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lie in over-dense regions. 

Several approaches have been used to investigate the 
nature of high— z DLAs. Broad band imaging surveys of 
DLAs have been used since Steidel & Dickinson (1992); 
Steidel & Hamilton (1992) were unsuccessful in finding 
the absorber, but did find large numbers of Lyman break 
galaxies (LBGs) (e.g. Steidel & Hamilton 1993). In addi- 
tion, it has become clear that DLAs (1) show little evo- 
lution from z = 4 to z = 1 in their contribution to the 
matter density of the Universe fljnhAiz) (e.g. Rao & Turn- 
shek 2000; Boissier, Peroux & Pettini 2003, and references 
therein); (2) show little metallicity evolution (Pettini et 
al. 1999; Prochaska & Wolfe 2002) from z = 3 to z = 0.3; 
(3) have heterogeneous chemical properties (Pettini et al. 
2000); (4) have various morphologies (Le Brun et al. 1997; 
Rao & Turnshek 2000); and (5) have low star formation 
rates (Kulkarni et al. 2000; Bouche et al. 2001). 

These results indicate that DLAs do not participate sig- 
nificantly in the overall chemical enrichment of the Uni- 
verse and hence do not trace star formation in any straight- 
forward way. Both these observations and theoretical work 
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show that DLAs are not a single type of galaxy, but are 
merely characteristic of a type of region: namely, one with 
a large neutral column density (e.g. Khersonsky & Turn- 
shek 1996). 

Similarly to broad-band imaging, early attempts to de- 
tect diffuse Lya emission from DLAs at 2: > 2 using 
deep narrow band imaging (e.g. Lowenthal et al. 1995) 
did not reveal the absorber, but unveiled a few companion 
Lya emitters, hinting at the clustering of galaxies around 
DLAs. This prompted Wolfe (1993) to calculate the two- 
point correlation function at z = 2.6 and to conclude that 
indeed Lya emitters are clustered near DLAs. Some re- 
cent Lya searches succeed in unveiling the absorber (e.g. 
Fynbo et al. 1999). 

Here, we seek to constrain the matter distribution 
around DLAs at 2 ~ 3 using LBGs as tracers of the 
large scale structure. Specifically, are the DLAs in over- 
or under- dense regions? Previous attempts at answering 
this question include the investigation of Gawiser et al. 
(2001), which ruled out any clustering between galaxies 
and a DLA at 2; ~ 4. Recently, Adelberger et al. (2003) 
found only two galaxies within 265"(~ 5.7h~^ Mpc) and 
within |Az| < 0.0125 (- 8^."^ Mpc) of four DLAs at z - 3 
whereas 5.96 were expected. They conclude from the same 
null hypothesis test used here (see below) that the lack of 
galaxies near the four DLAs is significant at > 90% level 
and argue that this is evidence that DLAs and LBGs 'do 
not reside in the same parts of the universe'. 

In contrast, Francis & Hewett (1993) reported the dis- 
covery of super-clustering of sub-DLAs at 2; ~ 2.3. Recent 
narrow-band imaging (Fynbo et al. 2003, private com- 
munication) of the same field shows that spectroscopically 
confirmed Lya-emitting galaxies are clustered at the red- 
shift of the two strongest Hi clouds. These observations 
and others (e.g. d'Odorico, Petitjean & Cristiani 2002) 
show that at least some DLAs reside in overdense regions. 

The field of QSO 3C336 shows an over-density of galax- 
ies around an intermediate redshift DLA at z — 0.656 
(Steidel & Dickinson 1992; Steidel et al. 1997; Bouche et 
al. 2001). Steidel & Dickinson (1992) suggested the pos- 
sibility of a cluster at z = 0.656. From a sample of spec- 
troscopically confirmed galaxies within 50" (0.4 Mpc 
co-moving) of the QSO line of sight, Steidel et al. (1997) 
showed that galaxies are clustered in redshift around the 
metal-line systems (see Steidel et al. 1997, Fig 7). A closer 
analysis reveals that they are even more clustered around 
the DLA: the number of galaxies per unit redshift, dN/dz, 
peaks at z = 0.656 (aside from the peak at Zem = 0.927). 
Indeed, out of 11 galaxies at 0.4 < z < 0.85, four galaxies 
lie within < 0.03 (9,000 km s^^) of the DLA on scales 
of ^ 200/i~^ Mpc (co-moving). 

In this paper, we present results of our survey of 
LBGs around a DLA at ^ = 2.974 towards QSO APM 
08279+5255. This field is the first of four. The results 
on the other three (one with some spectroscopy) are still 
pending and will be presented in Bouche & Lowenthal 
(2003) along with the details of our procedure and data 
analysis. Our study is similar to the investigations of Ga- 
wiser et al. (2001) and Adelberger et al. (2003); however, 
our MOSAIC sky coverage is more than 25 times larger, 
providing us with both extended spatial sensitivity and a 

^ Their study was aimed at investigating Lya forest between 3.109 < 



built-in control sample of the background surface density 
of LBGs. 

Throughout this paper, we adopt = 0.3, Qa = 0.7 
and Ho = lOO/i km s~^ Mpc~^; thus, at z ~ 3, 1" corre- 
sponds to - 21.5/i"i kpc and 1' to ^ 1.29/1"^ Mpc, both 
co-moving. At that redshift, H{z) ~ AAQHo, so 6z = 0.1 
corresponds to 67h~^ Mpc in co-moving coordinates. Mag- 
nitudes are Vega based unless stated otherwise. 

2. THE DATA 

2.1. The APM 08279+5255 field 

This field, like the other three in our survey, was selected 
for the presence of a DLA at 2; ~ 3 and with requirement 
Zabs *C Zem- The rcdshift range z 3 is ideal for selecting 
LBGs efficiently using standard photometric filters. 

Even taking into account its gravitational lens magnifi- 
cation 10 - lOOx, this APM 08279-^5255 ranks among 
the most luminous objects in the Universe (e.g. Irwin 
et al. 1998; Lewis, Robb & Ibata 1999). The DLA at 
Zabs = 2.974 was first reported by Petitjean et al. (2000) 
based on a high resolution spectrum (6 km s~^) of the 
QSO APM 08279+5255 (zom = 3.91), obtained by EUison 
et al. (1999) ^ at the 10-meter Keck telescope with HIRES. 
Due to uncertainty in the continuum, the DLA's column 
density is poorly constrained: 19.8 < logA'"(Hl) < 20.3. 
The absorber could therefore technically be a 'sub-DLA'. 
However, the Lya equivalent width VFrcst(Lya) > 4.8A 
(Petitjean et al. 2000) and its strong Mgii equivalent width 
Wrest (Mgii) > O.eA (Kobayashi et al. 2002) supports its 
classification damped cloud. 

Petitjean et al. (2000) estimated a minimum cloud size 
tdla > 200/175^ pc from the image separation of the QSO. 
The most robust estimate is from Tzanavaris & Carswell 
(2003) who found Tdla > 450/i~^ pc (assuming the grav- 
itational lens is at z 1). The QSO and DLA properties 
are summarized in Table 1. 

2.2. Observations, data reduction & completeness 

The observations, summarized in Table 2, were carried 
out with the MOSAIC camera (Jacoby et al. 1998) at 

the Kitt Peak National Observatory 4-m telescope on UT 
February 7th and 8th, 2000. The wide field imager MO- 
SAIC has eight 2kx4k thinned SITe CCDs. With 0.258" 
per pixel, it has a field of view of 36 arcmin on a side. 

Our field was imaged through the UBV&iI filters with 
integration times as shown in Table 2. The four nights 
were photometric. A detailed description of the data re- 
duction and calibration of our survey is deferred to Bouche 
& Lowenthal (2003). To summarize, the images were bias- 
subtracted, corrected for the MOSAIC ghost pupil image, 
dome- and sky-fiatfielded, cleaned for cosmic-rays, depro- 
jected and registered onto a common frame, before being 
CO- added. All magnitudes quoted below are measured in 
a 2 X FWHM diameter aperture, where the FWHM is 
the seeing. The resulting magnitude limits are shown in 
Table 2. 

From Monte-Carlo simulations, we estimated our com- 
pleteness as a function of magnitude, and is shown in 
the inset of Fig. 1 as the dotted line. Our 50% com- 
pleteness level is m/ = 24.35. Using the transformation 

< 3.701. 
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to convert our /-band magnitudes into AB magnitudes, 
i.e. Iab = iTii + 0.47, tliis corresponds to Iab = 24.8 
{TZ ~ Gv) and to ~ 0.6L*, where we used = 24.5 
at ^ ~ 3 (Steidel et al. 1999). In terms of L* (today), 
it corresponds to L/L* ~ 3 using a distance modulus 
m — M = 46.25 — 51og/i at z = 3 in our cosmology, 
M* AB -5 log/i = -20.04 from SDSS (Blanton et al. 2001), 
and a fc— correction between our / and the SDSS g filter of 
—0.30 (assuming the Irr SED of Coleman, Wu & Weedman 
1980). 

2.3. Sample selection 

Sources were detected in the J-band using SExtractor 
(Bertin & Arnouts 1996), and aperture magnitudes were 
measured in all four bands using identical matched aper- 
tures defined in the /-band image to extract 30,000 objects 
with / > 22.5. Fig. 1 shows the position of each object in 
the {U — B), (B — I) color space. We first selected a subset 
(1/3) of our /-band catalog using the color cuts shown by 
the dot-dashed lines in Fig. 1. This removes most of the 
z < 1 objects and reduces the amount of computing in 
the next step. We then estimated photometric redshifts 
of all objects in our subsample using the template fitting 
routine Hyperzbom Bolzonella, Miralles & Pello (2000) 
This constrains the redshift mainly via the shape of the 
SED at Xrest = 912A. In addition, we used the appar- 
ent magnitude as a redshift prior following the Bayesian 
prescription of Benitez (2000). The advantage of this last 
step is that it breaks most of the degenerancies of the 
SED fitting procedure (i.e. 4000A break taken as a Lya 
break) as shown by Benitez (2000). Finally, from our cat- 
alog now complemented with photometric redshifts ^;phot 
and likelihood distributions P{z), we selected galaxies with 
22.5 < I < 24.35 and computed the probability of each to 
be in three redshift slices of width Wz = 0.15: (1) Pdla = 
P(^DLA ± 0.075), (2) P+ = P([zDLA + 0.15] ± 0.075), (3) 
P_ =P([^DLA- 0.15] ±0.075). 

Out of 428 LBG candidates at redshift 2.75 <z < 3.25, 
84 objects with Pqla > 0.5 constitute our sample of can- 
didates in a redshift slice centered on the DLA and are 
shown in Fig. 1 as open squares. This 0.5 threshold maxi- 
mizes the number of objects and minimizes the likely out- 
liers: PoutUer = Pz[\^z\ > 0.2(1 + z)] < 0.1, foUowiug 

Benitez (2000). As a posteriori check, the distribution of 
LBG candidates in Fig. 1 shows that we are not excluding 
significant 2 ~ 3 objects in our color pre-selection. The 
number counts arc shown in the inset of Fig. 1. Similarly 
89 (17) LBG have P+ > 0.5 (P_ > 0.5), respectively. 

The slice width {Wz = 0.15) was chosen to correspond to 
the rms of the residuals between photometric and spectro- 
scopic redshifts ^■(A^). ^(A^) is ~ 0.15 ((t(Az/(1 + z)) ~ 
0.06) in the redshift range 2.8 < z < 3.5 from our analysis 
of photometric redshifts in the HDF-N. The residuals are 
similar to the accuracy in other studies (e.g. Fernandez- 
Soto et al. 2001). Wc defer the detailed description of 
our photometric redshift analysis to Bouche & Lowenthal 
(2003). 

The slice width adopted here produces the largest sam- 
ple in the smallest redshift slice. Our results are not af- 
fected when we chose a slice of half or twice the current 
size. At 2; = 3, the width Wz = 0.15 corresponds to 

* We note that we used the MOSAIC filter curves and the CCD resi 



~ 100 Mpc (co- moving) or 25 Mpc (physical) assuming 
h = 0.75. Pecuhar velocities of 7,500km would be re- 
quired to have an effect on our results. This figure is much 
larger than any peculiar velocity in local clusters of galax- 
ies. At z = 3, peculiar velocities are likely to be smaller 
than 1,000km s~^. Therefore, peculiar velocities are not 
affecting our results. 

3. RESULTS 

3.1. Clustering analysis 

Hierarchical models of galaxy formation predict that 

massive galaxies will likely be found in regions of high 
density, whereas low-mass galaxies are more uniformly dis- 
tributed. This produces an enhancement of the clustering 
of high-mass galaxies since they tend to form near each 
other in rare overdense regions. Therefore, clustering is a 
probe of the mass distribution of galaxies. For instance, if 
LBGs cluster around DLAs, it would indicate that DLAs 
reside in massive dark halos. 

A natural tool to investigate the clustering of galaxies is 
the two-point auto-correlation function £,gg{r) = {r/ro)~^. 
However, our redshift accuracy ^(A^) ~ 0.15 corresponds 
to lOO/i^^ Mpc. Therefore, without more accurate red- 
shifts, the angular auto-correlation function Wgg{0) is a 
more relevant statistic. Similarly to the two-point angu- 
lar correlation function Wgg , the cross-correlation between 
DLAs and LBGs, cudg, gives the excess probability of find- 
ing an LBG in a small annulus of area dn between 9 and 
+ d6 from a known DLA. The expected number of LBGs 
in that annulus is then: 

N{LBG\DLA) = £,(1 + LVdg{0))dn, (1) 

where Eg is the imderlying surface density of LBGs which 
is = 0.041/i^ Mpc^^ in our slice centered on the DLA. 
Thanks to the large field of view of the MOSAIC imager, 
there is little ambiguity in estimating the background sur- 
face density of LBGs. The background density is measured 
at 15 < r < 20h~^ Mpc, i.e. where galaxies are physically 
unrelated to the DLA. 

Estimating w^g accurately would require a large number 
of DLAs. Therefore, we investigated the spatial distribu- 
tion of LBG candidates in two limiting cases: (i) assuming 
the LBG-DLA cross-correlation is the same as the LBG 
auto-correlation and (ii) assuming the cross-correlation is 
zero. In addition, we used a 3-D estimator of the cross- 
correlation. In the next two sections, we present our re- 
sults and discussion of this analysis. 

3.2. Clustering results 

Figure 2 (left panels) shows the spatial distribution of 
LBG candidates in the three redshift slices, each of width 
Wz = 0.15: a) the slice at z = zdla — 0.15 with P_ > 0.5, 
b) the slice centered on the DLA at zdla = 2.974, and c) 
the slice at z = zdla + 0.15 with P+ > 0.5. We smoothed 
the surface density with a Gaussian kernel of 2.5/i~^ Mpc, 
shown with the contours in Fig. 2. The contours show the 
surface density S = (S— < S >)/ < S >, where < E > is 
the average surface density of the field. The thick contour 
shows 5 = 0, the continuous contours show 6 = 1, 2, 3, and 
the dotted contour shows 6 = —1. 

curves in performing the SED fits. 



4 



The right panels of Fig. 2 show the radial cumulative 
surface density centered on the DLA for the same three 
slices. The error bars are derived using bootstrap resam- 
pling method Mo, Jing & Boerner (1992). 

The middle panels of Fig. 2 show that the surface den- 
sity of LBG candidates is higher by a factor ~ 3 x within 
2.5 — 5 Mpc than the averaged background galaxy density 
T,g. This shows that at least some DLAs could reside in 
overdense regions, and implies a possible cross-correlation 
between this DLA and LBGs. 

For the slice centered on the DLA, Fig. 2(b), we used 
the estimator of the cross-correlation proposed by Eisen- 
stcin (2002), which gives the volume average of the cross- 
correlation < >. This estimator uses weightings that 
produce a spherically symmetric aperture rather than the 
cylindrical one used above. Fig. 3 shows < > as a 
function of the scale length of the Gaussian window. We 
find a 3-D overdensity 1+ < >= n/rig ~ 4.97(±2.99) 
using a Gaussian window of 5h~^ Mpc. The errors include 
the shot noise and the contribution of the clustering of the 
galaxies. We again conclude that this DLA seems to reside 
in an overdense region. 

Next, in order to assess the significance of our result, we 
examined the two limiting situations presented and some 
possible sources of contamination. 

In the first case, where wc assume the LBG-DLA cross- 
correlation is the same as the LBG auto-correlation, i.e. 
^dg = ^gg, WC tcst whether this excess of galaxies around 
the absorber could be caused by the LBG clustering. A 
qualitative way to answer this is given by the dotted line 
in Fig 2, which shows the expected number of galaxies 
around the DLA due to the clustering of LBGs. That 
number is naturally found from Eq. 1 integrated from to 
^max ^jgj^g f^iyf, known auto-correlation of LBGs, £_gg, from 
Adelberger et al. (2003) converted into U!{9) using our red- 
shift distribution. Clearly the angular clustering of LBGs 
alone cannot explain our data. 

A more quantitative answer is given by computing the 
expected number of galaxies iVcxp('") again assuming S^dg = 
^gg , and comparing it with the observed number of galaxies 
iVobs- Table 3 shows the values for at various radii. E.g., 
within 2.5 < re < 5h~^ Mpc, we find 7 galaxies while 2.70 
were expected. Could Nohs and A/exp be drawn from the 
same Poisson distribution of true mean iVcxp? The answer 
is given by testing the null hypothesis "iVobs and iVexp are 
drawn from the same Poisson distribution" versus 'Wobs is 
greater than A^cxp" , i-e. the observed overdensity is statis- 
tically significant. At 2.5 < rg < 5h^^ Mpc, we find there 
is a probability of finding Nobs > 7 galaxies (under the null 
hypothesis with Afexp= 2.70) of < 1% (=P-value shown in 
Table 3). This implies the null hypothesis ought to be re- 
jected, and the observed number of galaxies is larger than 
what one would expect from the LBG auto-correlation at 
the > 99% significance level. 

In the second limiting case Udg = 0, we test whether 
a purely random fluctuation in S could cause our re- 
sult. This would decrease the expected number of galaxies 
slightly and therefore increase the significance level (see 
Table 3). 

3.3. Discussion 



Our results are not sensitive to the background den- 
sity. The background density, itself, is a strong function 
of the magnitude limit and of the redshift slice width. A 
change of the background density will change signifi- 
cantly both the observed and expected number of galaxies 
throughout the field. However, both are affected in the 
same way from Eq. 1 by Eg. The effect on the overdensity 
cancels out, since the overdensity is the excess of observed 
to the expected number of galaxies. The same relative 
overdensity by a factor of ^ 3x is observed with different 
magnitude cuts from 24 to 24.5, and with a width half or 
twice our current redshift width. 

Alternative causes for the apparent excess of galaxies 
include contamination from low redshift sources (stars or 
galaxies) and signal from the cross-correlation between 
LBGs and other quasar absorption line clouds besides the 
DLA. 

We estimate our contamination from low-redshift galax- 
ies to be ~ 10% at 2.8 < z < 3.5 based on our photometric 
redshift analysis in the HDF (Bouche & Lowenthal 2003). 
Another source of contamination is from faint stars with 
color similar to those of our LBG candidates. To estimate 
the stellar contamination, we used our star coimts up to 
mi < 23.5 and extrapolate them to our magnitude limit, 
assuming the counts to be constant at m/ > 23.5. This 
gives an upper limit on our stellar contamination of < 7%. 
We note that both of these contaminations would dilute 
any clustering signal. 

A high-resolution spectrum of this QSO (Ellison et al. 
1999) shows there are 12 Civ metal-line systems with 
3.1< z <3.5 and 4 with 2.97< z <3.10 (Tzanavaris 2002). 
However, the cross-correlation between Civ systems and 
LBGs is less strong {rg = 3h~^ Mpc) than the LBG 
auto-correlation (fo = 4/i~^ Mpc, Adelberger et al. 2003). 
Therefore, on scales 5h~^ Mpc, the cross-correlation be- 
tween Civ systems and LBGs cannot bias our results, for 
it would have to be larger than the LBG auto-correlation. 

We close this section with a discussion of our search for 
the DLA itself. The QSO is not only lensed, but it also 
saturated the detector, thus preventing us from any at- 
tempt at PSF subtraction. Due to the brightness of the 
QSO, we are not sensitive to faint objects with impact 
parameters less than 5-6". From our source catalog, we 
searched for all objects within 45" (285 kpc in physical 
units) that have colors, i.e. photometric redshift, consis- 
tent with being at z 3. Wc find one object at an impact 
parameter of re = 145 kpc (physical) that did not make it 
into our final catalogs for it is fainter than our complete- 
ness limit by 0.3mag. It has a magnitude of 24.65 ± 0.2 
(M/_AB = —21.35) and a redshift Zphot = 3.03. The colors 
of this galaxy are shown in Fig. 1 by the open triangle. Be- 
sides this galaxy, the closest object that did end up in our 
final catalog is 2.35 Mpc (co-moving) or 590 kpc (physical) 
away. We regard the galaxy 145 kpc away as an unlikely, 
but possible, candidate for the DLA. 

4. SUMMARY AND CONCLUSIONS 

To summarize, in the field of APM 08279+5255, we 
find 7 galaxies with re < 4/i~^ Mpc (186") and within 
jA^I = 0.075, whereas one would expect 2.70 assum- 
ing that the cross-correlation is given by the LBG auto- 
correlation, i.e. that DLAs would be galaxies similar to 
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LBGs. The hypothesis that this is due to either the galaxy 
autocorrelation or to random fluctuation is rejected at the 
> 95% level. We estimated our contamination from stars 
or I0W-2; galaxies to be < 17%, since such contamination 
would only dilute our cross-correlation, so our results gives 
a lower limit on the cross-correlation. 

We therefore conclude that^ compared to the averaged 
background galaxy density T,g, there is an overdensity 
S/Sg ~ 3 of LBGs within 2.5 — 5 Mpc around this 
2aiw = 2.974 DLA. This scale corresponds to only 1.3 Mpc 
in physical coordinates at z = 3, when clusters would still 
be in formation and have core radii larger than 8 Mpc. 
Large scale structures at high redshifts shouldn't be un- 
expected. Recently, Shimasaku et al. (2003) reported the 
discovery of large-scale structure of Lyman alpha emitters 
at 2; = 4.86, 12h^Q Mpc (co-moving) in radius, based on 
wide-field imaging with Suprime-Cam on the Subaru tele- 
scope. 

Our result implies a possible correlation between this 
DLA and LBGs, and that at least some DLAs reside in 
overdense regions. However, preliminary results in two 
of our other fields (Bouche & Lowenthal 2003) indicate a 
weak and an anti-correlation between LBGs and the other 
two DLAs. This could imply that not all DLAs reside in 
the same environment. This reinforces the idea that DLAs 
do not likely represent a single type of object, but are a 



wide variety of objects with high neutral hydrogen column 
density. 

As a final note, the galaxy with the smallest impact 
parameter in our final slice centered on the DLA is 
2.35hr^ Mpc (co-moving) or 590 kpc (physical) away from 
the QSO line of sight. In addition, we search for all objects 
within 45" (285 kpc physical) and find an LBG candidate 
with Zphot = 3.03 at = 145 kpc, with a luminosity of 
Mj^AB — —21.35. Due to the large impact parameter, this 
galaxy is unlikely to be responsible for the DLA, although 
it may be associated with it. 

We are currently obtaining multi-object spectra with 
Gemini+GMOS of our LBG candidates within 5 Mpc 
of the QSO line-of-sight to investigate further the cluster- 
ing of LBGs and the DLA. 
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Table 1 
Properties of the field. 



Parameter Value 



QSO APM 08279+5255 {z^m = 3.91) 



R.A. (J2000) 
Decl. (J2000) 
niR (mag) . . 
^c/'*(mag) .. 

Damped Ly-a Cloud (^;abs = 2.974) 



logNHi^ 19.8-20.3 

[Fr/Ilf' -2.:-!l 



08''31™41.6*' 

52°45'17" 

15.2 

0.20 

0.04 



References. (a) from Schlegel et 
al. (1998), averaged over the field; (b) 
Petitjean et al. (2000) 



Table 2 



Summary of observations. 





Exp. /Frames 


Airmass 


FWHM 


/i/(l(T) 


m7(3a) ^ 


Completeness 50% ^ 


Filter 


(sec./#) 


(min-max) 


(arcsec) 


(niAB /arcsec^) 


(hiab) 


(mAs) 


U 


13500/15 


1.07-1.17 


1.1 


28.61 


26.68 




B 


2100/7 


1.14-1.27 


1.1 


28.37 


26.42 




V 


3000/10 


1.21-1.37 


1.2 


28.23 


26.19 




I 


7590/20 


1.11-1.62 


1.1 


27.62 


25.68 


24.80 



measured inside a 2 x FWHM diameter aperture. 



Table 3 

Null Hypothesis Test at various radii. 



re {h~^ Mpc) 


S {h^ Mpc-2) 




N ^ 


P-value'^ 


0-2.5 


0.05 


1 


1.04 


0.27 (0.19) 


2.5-5 


0.12 


7 


2.70 


<0.01 (<0.01) 


5-10 


0.07 


15 


9.90 


0.10 (0.07) 


15-20 


0.04 


16 


16.71 


0.50 (0.43) 



^Observed number of neighbors. 

''Expected number of neighbors assuming = ^gg. 

'^P-values of the null hypothesis, i.e. probability that Nohs>N^ 
if iVobs were drawn from a distribution of true mean equal to N^. 
assuming ^dg = ^gg (assuming ujdg = 0). 
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Fig. 1. — Color-color {U — B) ab , {B—I)ab diagram. Each dot represents one of our ~ 30, 000 objects in our image with 23.0 < Iab < 24.80. 
Open squares are our LBGs candidates in the redshift slice centered on the DLA {z^bs = 2.974), i.e. with Pdla > 0.5 (see text). Objects 
with U and B upper limits (filled triangles) have their U — B set to an arbitrary constant (= 5), and the B magnitude limit is used in S — /. 
Stars show the locus of stellar objects with 18.5 < Iab < 20.5. The open triangle is the closest LBG to the line of sight, which has an impact 
parameter of 145 kpc (see § discussion). The evolutionary track of an Irr SED (averaged over different dust content, from Av = to Ay = 1.2) 
is shown in redshift steps of 0.1 for illustrative purposes; the z = 3 mark is shown with the large filled square. We pre-selected objects within 
the color cut shown with the dot-dashed lines. The inset shows their number counts (AT) as a function of magnitude (Iab)- The dotted line 
shows our completeness (/). 
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APM 08279+5255 

(a) 
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Fig. 2. — The spatial distribution of LBG candidates (left) and the radial surface distribution (right) for three different redshift slices: a) 
slice centered at z = zola + 0.15, with P+ > 0.5; b) slice centered at 2 = zdlAi with Pdla > 0.5; and c) slice centered at z = zdla — 0.15, 
with P— > 0.5. Panel b) is centered on the DLA, marked by the cross. Panels a) and c) shows the distribution of galaxies unrelated to the 
DLA. Contours on the left panels shows the surface densities <5 = S/S — 1, where S is the mean density over the entire field. The thick 
contour shows <5 = 0; the continous contours show 5 = 1,2,3; the dotted countours show 5 = —I. The error bars on the right panels are 
computed using bootstrap resampling. The QSO line-of-sight is maked with a cross. For galaxies within the slice centered on the DLA (b), 
the surface density distribution shows a peak (significant at > 95%; see text) at ~ 4 Mpc from the DLA. This indicates an over-density of 
galaxies by a factor ~ 3x on scales 2.5 — 5 Mpc from the DLA compared to the averaged background galaxy density Eg. 
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Fig. 3. — The volume averaged (using a Gaussian window function) of the cross-correlation function ^^g > following Eisenstein (2002) is 
shown as a function of the scale length a of the window function. The (Icr) errors, shown as the dotted lines, include the Poission variance 
and the variance from the clustering of the galaxies. 



